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Abstract — We consider a list decoding algorithm recently 
proposed by Pellikaan-Wu [8] for g-ary Reed-Muller codes 
■Si^£q(£, m, n) of length n < q'" when £ < q. A simple and 
easily accessible correctness proof is given which shows that 
this algorithm achieves a relative error-correction radius of 
t < y 1 — \J lq m ~ x /nj . This is an improvement over the proof 
using one-point Algebraic-Geometric codes given in [8]. The 
described algorithm can be adapted to decode Product-Reed- 
Solomon codes. 

We then propose a new low complexity recursive algebraic 
decoding algorithm for Reed-Muller and Product-Reed-Solomon 
codes. Our algorithm achieves a relative error correction radius 

of t < 175=1 f 1 — ^ki/qj. This technique is then proved to 
outperform the Pellikaan-Wu method in both complexity and 
error correction radius over a wide range of code rates. 

I. Introduction 

With the discovery of deterministic list-decoding algo- 
rithms for several Algebraic-Geometric codes, most notably 
the Guruswami-Sudan [6] algorithm, there has been renewed 
interest in algebraic decoding methods for other related q- 
ary codes such as the Reed-Muller [7], [8] and Product-Reed- 
Solomon [9] codes. However some of the existing correctness 
proofs for these algorithms use advanced algebraic geometric 
tools. In this paper we first derive a proof for a list decoding 
algorithm for a g-ary Reed-Muller code. Our proof is from 
first principles and require only the most basic notions from 
finite field theory. We then proceed to propose new recursive 
list decoding algorithms for Reed-Muller and Product-Reed- 
Solomon codes. These algorithms are rigorously shown to 
outperform the Pellikaan-Wu method in both complexity as 
well as error-correction-radius. 

The basic idea of our new proof for the Pellikaan- 
Wu algorithm is to "lift" a multivariate polynomial in 
F q [x\,X2, ...,x m ] to a univariate polynomial in F^IX] using 
a deterministic mapping rule. This in turn results in a higher 
total degree polynomial. The increase in degree will not be 
high enough to render our list decoding strategy for Reed- 
Muller codes useless at meaningful rates. A higher degree 
for the lifted polynomial means that this Reed-Muller code 
list decoding algorithm has a lower relative error-correction 
radius (as a function of the rate) than a comparable rate 
Reed-Solomon list decoder based on the Guruswami-Sudan 



algorithm. In the following section we describe the mapping 
rule and the decoding algorithm in some detail. 

In the final section we propose new algorithms for de- 
coding Reed-Muller and Product-Reed-Solomon codes. Our 
algorithm is more efficient than the Pellikaan-Wu method by 
approximately a quadratic factor. Furthermore it outperforms 
the Pellikaan-Wu algorithm in error-correction-radius over a 
wide range of code rates. 

II. Correctness of a List Decoding Algorithm 
Let us begin by defining a q-aiy Reed-Muller code. 

Definition 1 The q-ary Reed-Muller code ' q (£, m, n) of 

length n < q m is defined as the set of vectors given by: 

def 

@JZ q (l, m, n)={ [<p(Oi) <p(a 2 ) ••• <p(a„)] 

I <peW q [xi,x 2 , ...,x m ], deg(<p) <l } (1) 

where {ffi, «2, • • ■ , <*„} are any set of n distinct points in F™. 
Here by deg(<p) we mean the total degree of the multivariate 
polynomial (p. 

The following well known property will be useful: 

Proposition 1 Let {a\, 02, ■ ■ ■ , a m } be a basis for ¥ q m over 
Fq and let [x\X2 ■ ■ ■ x m ] £ F™. Then the map \j/ : F^ —> V q m 
defined as in ([2]) is an isomorphism. 

def m 

[x\ x 2 ... x m ] h-> X = Y< a j x i ( 2 ) 

;=i 

For example one might as usual use a polynomial basis 
{l,^,^ 2 ,..., } where t, is any primitive element in V q m 



or even a normal basis of the form {£, C, q , ... 
where £ is a suitable primitive element in ¥ q m. 
Therefore we arrive at this elementary conclusion: 



C 9 }, 



Lemma 1 Let X £ ¥ q m. The reverse isomorphism for ([2]) is: 



X 

where 



[XIX2 



def 
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(4) 



is a non-singular (invertible) square matrix. 



OUTPUT: _§? 



Proof: Since X — Y!j=i a j x j, an d x j G ^q> we g et = 
YIj=\ a j x j using Fermat's little theorem. It only remains to 
show that A is non-singular. Note that in general A is not 
a Vandermonde matrix. However by construction, the set A. Complexity of Algorithm^ 
{oi, A2, • • • , fl m } is a basis for F^™ over F 9 . It then follows 
from [4, Corollary 2.38, pp. 58] that A is non-singular. ■ 

It follows from Lemma [TJ that there exist polynomials flj 6 
Fflm [X] of degree at most q'"~ 1 such that jc< = flj(X) ,l<j<m. 
Substituting for all xj in this manner, we have proved the 
following: 



This algorithm was originally proposed by Pellikaan-Wu in 
[8], though their proofs were different. 



The complexity of our proposed algorithm is of the same 
order as the complexity of Guruswami-Sudan algorithm for 
decoding Reed-Solomon codes over the extension field ¥ q m. 
This is ff(n^) field operations in ¥ q m. 



Theorem 1 Let n < q m . If £ < q then 

&Jg q {l, m, n) C My^n, V""')^ (5) 
where 3%y q m(n, £q m ~ l ) is the Reed-Solomon code given by 

^V(», V -1 ) = { [/(ft) /(ft) ... /(ft)] 

I /eF^degt/)^'"- 1 } (6) 

w/zere j3 ; = f £"' =1 ayG!^, ami a,- = f [a,i Ofe ■■■ C*fr?J> 1 < ' < 
n are f/ze points of evaluation for the Reed-Muller code. 
Moreover if the information polynomial associated with the 
Reed-Muller code is given by 

def m i 

<p(xi,X 2 , •••,*«) = £ < Pii,h,-MYl X 'j ^ 



i\ ,i2j— jim : 

Zjij<t 



7=1 



f/ien f/ie information polynomial f of degree at most lq m 1 
associated with the Reed-Solomon code is: 

m 

f{X)= £ (/>/;,- n^'/CA-ii'- (8) 

'i .'2 '»' : ;'=i 



Let t/// (jc, represent the Hamming distance between the 
two vectors. Using Theorem [TJ and the Guruswami-Sudan 
algorithm [6] for list decoding a Reed-Solomon code, we 
have proved the correctness of the following deterministic list- 
decoding algorithm for Reed-Muller codes: 



Algorithm 1 (RM-List-1) 
INPUT: q,e<q,m,n<q m 
STEPS: 



r=[nr 2 ... r n ]e¥" q . 



\-^/tq m - l /n 



1. Compute the parameter t - 

2. Using Guruswami-Sudan algorithm find a list of code- 
words c G i%y q m{n, £q m ~ x ) such that dn(c, r) < t. 

3. For every c£if check if c gFJ : 

i. If NO then discard c from _£f. 

ii. If YES then check if c € MJ( q {i, m, n) : 

a. If NO then discard c from J£. 

b. If YES then keep c in the list ££. 

4. return 



B. Comparison to previous results 

The Pellikaan-Wu algorithm for decoding Reed-Muller 
codes by means of embedding into one-point Algebraic- 
Geometric codes was shown [8] to achieve an error correction 
radius of n (\ — \/l{q + l) m ~ l /n\ . It is interesting to note 
that the error-correction radius demonstrated herein is always 
larger than that suggested by the Pellikaan-Wu formalism 
employing Algebraic-Geometric codes. However we believe 
that the more important contribution of this paper is the readily 
accessible correctness proof which relies on just a few basic 
notions from Galois theory. 



def 



C. Product Reed-Solomon codes 

Product Reed-Solomon codes ^>^y qM (q n \ h,..., k m ) 
®? = i&S e ' q (q, kj) over ¥' q can be thought of as the set of 
vectors whose q 1 " coordinates consist of the q m evaluations 
of m-variate information polynomials with coefficients in 
¥ q and degree in the i th -variable Xi at most (k, — 1). m 
is usually called the dimension of the product code. Thus 
^My q .m{q m , h,...,k m ) is contained in ^Jl q i^ =x (k { - 
1), m, q m ). When 1) < q the list decoding algorithm 

given in Algorithm [T] may be used essentially without any 
modifications. Several Product-Reed-Solomon algebraic list 
decoders, including a similar method as sketched above are 
described in [9]. Using Algorithm [T] it is possible to achieve 
a relative error correction radius of (1 — \/Y!iL\ Pi), where 

def, , 

Pi = h/q. 



D. Zeros of Multivariate Polynomials 

From Theorem [T] it is clear that f(X) being of degree 
at most £q' n ~ l , has at most iq m ~ l zeros in ¥ q m, including 
multiplicities. Therefore a non-zero multivariate polynomial 
<p(jci, X2, ...,x m ) of total degree £ has at most £q m zeros in 
¥ q \ This gives the famous DeMillo-Lipton-Schwartz-Zippel[2] 
lemma for polynomials over finite fields. Note that the state- 
ment above appears to be stronger than the classical lemma 
in that this counts multiplicities too. Moreover the proof also 
appears to differ from the traditional expositions which use 
probabilistic arguments. 

Next we propose a lower complexity recursive algebraic de- 
coder which outperforms the Reed-Muller decoder considered 
in this section. 
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Fig. 1: 2D PRS code described on a rectangular array. 



III. A Recursive Decoding Algorithm for 
Reed-Muller and Product Reed-Solomon Codes 

def 

For simplicity, let n — q m . A codeword in the code 
^SiS^qfn (q m , ki,...,k m ) can be described within an m- 
dimensional cube of side length q. See Figure [T] Let a 
codeword c (correspondingly a received word, r) be so de- 
scribed. We will find it convenient to write this vector as 
(Cj-j h i m }, where each of the indices i; take values in the 
range {l,...,q}. We further use the notation (c^j- J+ j'."', 1 *) 
to denote the (j — 1) -dimensional vector formed out of 



(c,-j j, ,,„]> when the coordinates indexed by (£;,!; 
are fixed at (aj,aj 



+ 1 j ■ • • i *'m) 

,a m ) and the rest of the indices are 
free. By the nature of the product code, {c^'^ + •'. "' "' } belongs 

to .S t'./ H l ,(,,'' '. A, kj .;. 

Now consider the following decoding algorithm for the code 

£P S$,5^ q^mitf ) ^1 3 • • • ) Kn ) - 



def 



Algorithm 2 (PRS-Decoder) 

INPUT: q,(k\ 7 k2,. . . ,k,„) : kj < q,m; r <E W q , where r 

(K,/ 2 ,...,iJ);l <ij<q- 
STEPS: 

1. If m—\ do: 

i. Compute the parameter t\= q { \ — \Jk\jq^ . 

ii. Using Guruswami-Sudan algorithm find a list S£\ of 
codewords C\ € ^%S fi q {q 1 k\) such that dn{c\, { }) < 
h. 

iii. Search S£\ for C\ such that dnipx, (r,-,}) is least. 
Substitute in-place the positions corresponding to (r;. } 
in r with C\ and return . 

2. For a m = 1,2,..., q do: 

i. Set r 1 <-(r? m ; ,■ } 

1 <1,«2,— ,'m-l y 



ii. Sef m' <— m—\ and n' *— q m 

iii. Recursively decode r 1 using PRS-Decoder with input 
parameters q,(k\,k2,-.., k m > ) , m' ; r' G . 

3. Compute the parameter t m = q ^1 — \Jk m j q^ 

4. For each m— I tuple (ai,a2,...,a m -i) do: 

i. Using Guruswami-Sudan algorithm find a list 
J£ m of codewords c m £ &S^ q {q, k m ) such that 
d H (c m ,(r^' a2 - a "-'}<t m . 

ii. Search J£ m for c m such that dn{c m , |/-^ ,a2 '""" m - 1 1) 
/east Substitute in-place the positions corresponding 
to{r^ a2 '-'"'-'}withc m . 

5. return 

OUTPUT: Resulting vector r 

The following recursive algorithm uses PRS-Decoder to 
decode Sfi^tt q {Jt, m, n). 

Algorithm 3 (RM-List-2) 

INPUT: q,£ < q,m,n < q m ; r = [n r 2 . . . r n \ € ¥\. 
STEPS: 

1. For eac/z possible m-tuple (k\,hi,.. .,k m ) : kj < q,Y.jkj < ^ 
ao; 

i. Using PRS-Decoder with input parameters 
q,{k\,k2,-.., k m ) ,m; r £ F'J, decode r as c. 

ii. AaW c to a list Jzf o/ codeword candidates. 

2. return 
OUTPUT: «£f 

We have the following result concerning the decoding power 
of Algorithm [2] and Algorithm [3] 

Theorem 2 Algorithm^has a relative error correction radius 
of X m = f n"li (1 — \/Pi)' where p, = ki/q. Moreover, there exist 
error patterns of weight above n Il/!Li (1 — \fPi) which cannot 
be guaranteed to be efficiently decoded by Algorithm [3] 

Proof: 

Our proof is by induction. When m = 1, the claim is trivially 
true. Let us assume the claim to be true for some m — M. We 
will now show it to be true for the case m = M+ 1. Let there be 



a maximum of f^f+i = q n*=i (1 — -\/p7) errors. In Step 2 
of Algorithm |2j let there be a maximum of x recursions which 
fail to decode correctly. Since by the induction hypothesis, this 
would mean that there are more than tM errors in these x sub- 
recursions, we have that xtM < Substituting for f^+i and 
fM gives, x < q(l — ^/Pm+i)- These errors will get corrected in 
Step 4 of the algorithm. This proves the first part of the claim. 
For the 2D case, the proof is concisely depicted in Figure [2] 
To see the second part of the claim, we observe that 
an error pattern which is contiguously spread over an m 
dimensional sub-cube of volume more than «IT,li(l~ \/P~i) 
cannot be guaranteed to be efficiently decoded by the proposed 
algorithm. This shows that the error correction radius predicted 
in the first part of Theorem [2] is rather tight. ■ 



x rows are in error t\X = q ■ (1 — ^/p~i) ■ x 

*2 = <? 2 ■ (i - ypr) ■ (i - ypi) 

h<h =* * < 1 ■ (1 - Vft) 



a' 












































« 




a 

o 
































« 

* 




a 

o 

o 

a 











































Fig. 2: The proof of Theorem |5] for a 2D PRS code. 



A. Complexity of Algorithm |2] one/ Algorithm [3] 

Let # m be the complexity of decoding an m-dimensional 
Product-Reed-Solomon code using Algorithm [2] Then the 
complexity of decoding an m+ 1 dimensional code is # m +i = 
0(q&,„ + q'"$i). But #i = ^(^ 3 ) field operations in F 9 . 
This gives, #,„ = ff(q m+2 ) which is w ^(ra) for large m. The 
complexity of Algorithm b] is w &(n 2 ) field operations in F 9 . 
This is substantially better than the Pellikaan-Wu method in 
Algorithm [T] 

B. Comparison of Algorithm [7] and Algorithm [i] 

Algorithm [3] not only has a lower complexity, but also 
performs better over a wide range of rates. For example when 
£, p, > 1, the Pellikaan-Wu algorithm is not effective, whereas 
the new algorithm is still useful. Furthermore n?ii(l ~ y/Pi) 
is larger than (1 — Pi) for most code rates and the 

advantage is more pronounced at higher code rates. Figure 
[3] shows the decoding power of Algorithm [2] 

C. Other Related Product Code Decoders 

Several iterative hard decision decoders for Product-Reed- 
Solomon codes available in literature use some form of 
Algorithm [2] Usually such algorithms are described with 
no theoretical bounds on their error correction radii. These 
product code decoders find use in optical communication 
systems and LAN/WAN standards [10], [11]. Several hard- 
ware implementations of such decoders are commercially 
available[12], [13], [14]. Soft decision iterative decoders for 
product codes utilizing the "turbo-principle" have also been 




(a) Relative error correction radii of the two algorithms are 
compared when decoding a 2D Product-Reed-Solomon code. 
The surface which dominates for most of the rate region 
corresponds to the new algorithm. 




(b) 2D rate region where the new recursive 
algorithm performs better is shown. 
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(c) Relative rate region volume where our algorithm performs 
better is computed for various Product-Reed-Solomon codes of 
dimensionality m. The new algorithm is seen to out-perform the 
Pellikaan-Wu algorithm over much of the rate region. 

Fig. 3: Comparison of Pellikaan-Wu algorithm to our new 
recursive RM/PRS decoder. 



discussed in literature [5]. The performance of most of these 
hard decision iterative decoders can be very well characterized 
using Theorem [2] Similar conclusions are obvious for the case 
of other product codes which have algebraic bounded distance 
decoders available for their component codes. Theorem [2] 
implies the following for a general product code: 

Corollary 1 If for an m-dimensional product code P, there 
exists bounded distance decoders for each of its component 
codes such that the i th component code's decoder achieves a 
error correction radius oftj errors, then there exists a decoding 
algorithm for the entire product code P which can correct all 
errors up to a weight of t = YlT=i t i- 

The decoding algorithm for the code P mentioned in Corollary 
[T]can be obtained from Algorithm [2] with some obvious and 
minor changes and as such is not repeated here. This result is, 
to the best of the author's knowledge, the only such theoretical 
guarantee on the error correction radius of a general algebraic 
product code decoder. However for specific cases there are 
some stronger results available, for instance see the result of 
Lin-Weldon[l] for cyclic product codes. In another related 
example, Tanner[3] discusses bounds on a specific type of 
hard-decision decoder for product codes on graphs. In many 
cases of practical interest such bounds are difficult to apply 
because of their dependence on the knowledge of the girth of 
the underlying code graph. 

IV. Conclusions 

In this paper, we presented a simple and easily accessible 
proof for the Pellikaan-Wu algebraic list decoding algorithm 
for Reed-Muller codes. Our proof uses only the most funda- 
mental properties of finite field arithmetic. 



We also proposed a low complexity recursive algorithm 
for Reed-Muller and Product-Reed-Solomon codes. This new 
recursive algebraic decoding algorithm is then shown to have 
a significantly better error correction radius than the Pellikaan- 
Wu algorithm over a wide range of code rates. 
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